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Radio-frequency quantum engineering of spins is based on the dressing by a non resonant elec- 
tromagnetic field. Radio-frequency dressing occurs also for the motion of particles, electrons or 
ultracold atoms, within a periodic spatial potential. The dressing, producing a renormalisation 
and also a freeze of the system energy, is described by different approaches, dressed atom, mag- 
netic resonance semiclassical treatment, continued fraction solution of the Schrodinger equation. A 
comparison between those solutions points out that the semiclassical treatment, to be denoted as 
the S-solution, represents the most convenient tool to evaluate the tunneling renormalization of 
ultracold atoms. 

PACS numbers: 



I. INTRODUCTION 



The analysis of a system with few degrees of freedom, 
an electron or an atom, interacting with a large system, 
photons or phonons, relics often on a renormalization ap- 
proach, where the parameters of the initial system are 
modified by the interaction. Examples of this approach 
arc the effective mass for the electron motion in a semi- 
conductor, or the extensive renormalization in quantum 
electrodynamics. Another example is the dressed atom 
introduced by Cohen- Tannoudji in order to describe the 
modification of a two-level atomic magnetic response by 
an applied radio-frequency (rf ) field in the absence of de- 
coherence processes [l| . For a non- resonant rf driving at a 
high frequency Cohen- Tannoudji and Haroche [3| derived 
a renormalization of the atomic level splitting dependent 
on the amplitude of the rf field and described by a zero-th 
order ordinary Bessel function. That modification pro- 
ducing a magnetic "freezing" of the two- level response, 
(i.e. a nonmagnetic system), was examined for atoms 
in J3|-l6|, for a Bose- Einstein condensate of chromium 
in [J, for an artificial atom in Q , and recently proposed 
for improving the precision of optical clocks [9| . 

The same renormalization and freezing of the system 
properties under the application of a time-dependent 
modulation was applied to a variety of processes, all char- 
acterized by weak decoherence processes. We mention 
here the dynamical localization describing the renormal- 
ized motion of a charged particle within a periodic po- 
tential under a time modulated force [loj and the coher- 
ent destruction of tunneling for a double well potential 
with a periodic driving, with a complete localisation of 



a wave packet in one well for specific values of the driv- 
ing force [ll|. For the motion of ultracold degenerated 
atomic gases within a shaken optical lattice, the tunnel- 
ing atomic evolution is renormalized under the applica- 
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tion of a time modulated force, as proposed in jllMl 
and tested experimentally in jl7h23] within the frame- 
work of quantum simulation of solid state physics. The 
renormalization of an optical lattice potential acting on 
cold atoms in a regime of classical diffusion and transport 
was investigated in [23j . 

Target of the present work is to characterise the renor- 
malization and freezing for a wide parameter range. The 
starting point is the continued fraction solution for the 
two-level atomic magnetic response to an applied rf field, 
as typical of the nuclear magnetic resonance for a one- 
half spin. The regime of strong perturbation was investi- 
gated by several authors [24l - [28j . in presence or absence 
of decoherence. Their solution was expressed in terms 
of infinite continued fractions. The present work inves- 
tigates the renormalization process through the contin- 
ued fraction approach. That treatment allows us to ex- 
plore numerically the shaken-lattice renormalization for 
all parameter ranges, and in particular for explored ex- 
perimental conditions. The numerical complexity of the 
continued-fraction solution, and its slow convergence in 
the regime of experimental interest, brought us to con- 
sider carefully the corrections to the zero-order Bessel 
function derived for the dressed-atoni problem mainly 
through a semiclassical treatment |29l - [33 |. and later re- 
covered through renormalisation group techniques [33j . 
On the basis of the analogy between the renormaliza- 
tion of the magnetic resonance energy and of the atomic 
tunneling in optical lattice, we focus our analysis from 
the high frequency regime realized for a rf modulation 
at a very large frequency. Then we explore the correc- 
tions when this limiting condition is not precisely satis- 
fied. The low frequency regime covered by the treatments 
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of rcfs. [13, H, [3 nil is not examined here. 

Sec. II presents different systems where the renormal- 
ization of the interaction strength has been investigated: 
magnetic resonance, motion in a periodic potential, tun- 
nel coupling in a periodic potential. This section reports 
also the standard result of the renormalization process 
given by the zero-order Besscl function, valid under ap- 
propriate operating conditions. Sec. Ill reports the solu- 
tion for the temporal evolution of the wavefunction in the 
magnetic resonance case. Sec. IV derives the renormal- 
ization through the continued fraction approach, valid for 
all operating conditions, and also through a semiclassical 
treatment refening the Bessel-function result. Sec. V re- 
ports numerical results determining the limiting validity 
of the usual zero-order Bessel correction, and derive the 
renormalization for a large set of parameters. A conclu- 
sion completes our work. 



II. DRESSED SYSTEMS 

A. Magnetic resonance 

a) Semiclassical approach For a spin- 1/2 system inter- 
acting with a static magnetic field along the z axis and 
driven by an oscillating rf field along the x axis, the semi- 
classical Hamiltonian H^r is 



„ huJa 

tlsc — 



nn 



cos{iut)ax 



(1) 



where ax,z are the Pauli matrices and hujo the energy 
splitting between the magnetic levels and the Rabi 
frequency proportional to the rf field amplitude. 

By writing the atomic wavefunction \ip) written as a 
superposition of the |±) atomic cigenfunctions 



(2) 



The introduction of a renormalization coefficient R de- 
fined by the ratio between renormalizcd and original 
eigenvalues leads to 



(6) 



depending on the zero-order ordinary Bessel function 
JTq. Thus the applied sinusoidal magnetic interaction 
renormalizes the atomic coupling to the static magnetic 
field, with a reduction by the factor j7o(^). The effective 
magnetic energy is frozen whenever 51/a; is a zero root of 
the j/n Bessel function, as observed in the experiments of 
refs. 01 • The magnetic resonance renormalization 
was explored by ref. Q in the w < wq low-frequency 
regime, where the present approach is not valid. 

6 ) Quantized approach Introducing a quantum descrip- 
tion of the rf field, with operator and a for the cre- 
ation and annihilation of one radiofrequcncy photon, the 
dressed-atom Hamiltonian Hda of the above configura- 
tion is ^51 



where n represents the mean number of photons applied 
to the atoms. 

For the high-frequency case w ^ the last term 
in i/da niay be neglected and its eigenstates easily 
determined. Then a perturbation treatment for the Uz 
term of that Hamiltonian leads to E g. dp for describing 
the interaction with the static field [33. The previous 
renormalization result is obtained also through this 
approach. Notice that the dressed atom approach, 
and also the semiclassical approach of js^ , demonstrated 
that the j7o renormalization is valid for whatever spin 
value and equally spaced Zeeman levels. 



?n — ± 



the Schrodinger equation leads to the following temporal 
evolution for the Cm coefficients: 



iC+ =^C+ + ^cos{ijjt)C-, 



cos (ujt) C4 



(3) 



In the Wo — > limit, and for C+{t = 0) = 1 as initial 
condition, these equations have solution 



\m) 
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sin (cot) 
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(4) 

As in [7|,l27[ the result of calculating the az time- averaged 
mean value over jV') may be expressed through the fol- 
lowing renormalizcd eigenvalues of the Hsc Hamiltotian: 



Z Lu 



(5) 



B. Dynamic Localization 

Dynamic localization was introduced by Dunlap and 
Kenkre [ld| for the motion of an electron on a discrete 
one-dimensional periodic lattice with spacing in the 
presence of an oscillating force. It is based on exact cal- 
culations for the particle motion. A single-particle basis 
useful for describing the electron tunneling among the 
discrete lattice sites is provided by the j-th Wannicr func- 
tion centered on the j lattice site of the periodic poten- 
tial [3^ . In a given energy band the Hamiltonian for free 
motion on the periodic lattice is determined by tunnel- 
ing matrix elements, which in general connect arbitrarily 
spaced lattice sites. However, because the hopping am- 
plitude decreases rapidly with the distance, the tunnel- 
ing Hamiltonian may be well approximated by including 
only the hj tunneling energy hopping between neighbor- 
ing lattice sites. Under this hypothesis, the Hamiltonian 
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for the electron on the hnear lattice with an applied pe- 
riodic force Fcosiiot) is [lol [37| 

Hdi = hjY,i\]><3 + M + \] + ^><j\) 

m 

+ hKcos{Lot)Y,j\j ><j\. (8) 

3 

Here hK = Fdi^ is the time-modulated energy difference 
between neighboring lattice sites. Dynamic localization 
entails a suppression of the particle transport with the 
particle position oscillating in time and returning peri- 
odically to its original value. It is associated to particle 
motion on an infinite lattice and does not impose condi- 
tions on the frequency driving. Our focus based on the 
analogy with magnetic resonance is on the high frequency 
driving and on a lattice with a site finite number. 

For the case of two lattice sites (j = —1/2,1/2), in- 
troducing the Pauli operators, as = |l/2 >< 1/2| — 
I — 1/2 >< — 1/2| and so on, the above Hamiltonian be- 
comes f38| 

Hdi ^ hJa^ + ^cos{ujt)a^, (9) 

where we have introduced fl = K in order to emphasize 
the equivalence of this Hamiltonian with that of Eq. ([T]) 
apart a change of the quantization axes and the J = a;o/2 
parameter correspondence. Therefore the dressed atom 
renormalization applies also to this system, leading to a 
renormalized tunneling rate, 

Jeff = nj. (10) 

Once again the system response is frozen whenever 
il/uj = K/ihio) is a root of the Jq Bessel function. The 
applied sinusoidal force produces a dynamic localization 
of the particle. 

For a spin larger than one-half and more than two 
lattice sites, the Hamiltonian assumes a form equivalent 
to that of Eq. ([S]), except for the angular momentum. 
Thus the magnetic resonance analogy confirms the TZ 
renormalization also for an arbitrary number of lattice 
sites. 



C. Shaken optical lattice 

In a ID optical lattice ultracold atoms are confined 
within the potential minima created by a single laser 
standing wave with c?l spacing [39l . l40j . The Hamiltonian 
for atomic motion on the periodic lattice is determined 
by tunneling matrix typically including only the J hop- 
ping between neighboring lattice sites. 

A periodic force Fcos{ujt) (to be referred as lattice 
shake) drives the atoms inside the optical lattice. Us- 
ing the j-th Wannier function centered on the j lattice 
site, the Hamiltonian of Eq. ^ describes also the mo- 
tion of the ultracold atoms within the optical lattice. 



TABLE I: Values of the shaken-lattice experimental parame- 
ters, J equivalent to cjq/S , oj, Q = K/h, measured in units 
of the recoil frequency Urec/^Tv of the investigated atom. The 
last column reports the ujo/uj ratio. 



Ref. 


J = ioo/2 




S7 iJo /uj 




0.02-0.08 


0.15-0.9 


0-6 0.04-1.0 


[18j 


0.19 


0.8-4 


0-1.2 0.09-0.47 


[19j 


0.004 


1.9 


0-3 0.004 


[20] 


0.002 


0.9 


0-5 0.004 


[22] 


0.02 


0.9 


0-9 0.04 



with K ~ D, again the shaking energy difference between 
neighboring sites of the linear chain. Therefore the dy- 
namic localization and the renormalization of the previ- 
ous Subsection applies also to the ultracold atoms shaken 
lattices [TTI - [l6j , as tested in several experiments [T7l - {2^ . 
The parameters of the Hamiltonian of Eq. ^ investi- 
gated in those experiments are reported in Table I. Notice 
that most experiments investigated the high-frequency 
regime, but larg e deviations from that regime also oc- 
curred. Ref. [4l[ pointed out the difiiculties in the precise 
measurement of the tunneling freeze from the ultracold 
atoms images. 



III. CONTINUED FRACTION APPROACH 

If a single-particle Hamiltonian is periodic in time, 
with period T, then the Floquet's theorem [42| states 
the existence of a set of distinguished solutions \ipn{t)) 
to the time-dependent Schrodinger equation. These Flo- 
quet states, analogous to the usual energy eigenstates 
of time- independent Hamiltonian operators |2ll . [25l l43l | , 
have the form 

= M)) exp(-i£nt/n) . (11) 

with time periodic functions \u„{t)) ~ \un{t + T)). The 
quantum number n specifies the state. The quantities £„ 
are denoted quasicnergies. By inserting Eq. (|lip into the 
Schrodinger equation governed by the Hsc Hamiltonian, 
we deduce 

(^Hsc - ih-^^ \un{t)) = e„|u„(i)) , (12) 

to be regarded as an eigenvalue equation for the Floquet 
quasicnergies. The set of Floquet functions is complete in 
the Hilbert space on which acts the Hamiltonian. Hence, 
any solution \ip{t)) to the Schrodinger equation admits 
an expansion in the \un(t)) basis. 

If \un{t)) be a solution to the eigenvalue Eq. (fT^ with 
quasienergy e„, then |u„(<))e™'^* also is a T-periodic so- 
lution, with quasienergy e„ + mhuj m being an arbitrary 
integer, where ui = 2tt/T. Therefore the quasienergy of a 
Floquet state is determined only up to an integer multiple 
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of the hoj photon energy. In accordance with the sohd- 
state physics terminology, the quasienergy spectrum is 
said to consist of an infinite set of identical Brillouin 
zones of width fku, covering the entire energy axis, each 
state placing one of its quasienergies in each zone. 

The quasienergies may be determined by diagonaliza- 
tion of the Hamiltonian expressed in the Fourier space, 
as in [25| , or equivalently diagonalizing the dressed-atom 



Hamiltonian as in [35| . We will make use of the continued 
fraction solution of refs. [l^, H^- We apply the Fourier 
expansion to the C± coefficients of Eq. ^ 



n— +00 



The recurrence Eq. ([18)) has a continued fraction solu- 
tion pvl . [4^ , with expression for j > 



1 



(20) 



Dj + 2- — 



and similar expression for a negative j. By replacing xi 
and X-i into Eq. (|18p for j = 0, we obtain the continued 
fraction solutions for A+ and A_ , with 



A_ = 



1 



C-{t) 



iXt 



c+{t) 



iXt 



1 — — 00 1 



(13) 



Substituting these expansions into Eq. ([3]) and equating 
the same order Fourier components, one obtains 



A_ + ^ 1 + ^ 



il \ LJa j X 



and 



2 
2 



- nujj An = -— B„_i - — -B„+i, (14a) 
-nio] Bn = --An-i- -An+i. (14b) 



These equations can be separated into a first set with all 
even ^'s and odd B's being zero 



A, 



+ 



'Ak+i. 



(15a) 
(15b) 



and into a second one with all odd A's and even B's 
being zero 



4w 

(A--fc-)A. 



-Bk-i - B, 



k+l, 



(16a) 
(16b) 



All Floquet quasienergies are given by 



nhu} 



(22) 



(23) 



The continued fraction solution allows to determine nu- 
merically the Floquet quasienergies with the required ac- 
curacy. Fig. [T] reports the quasienergies within one Bril- 
louin zone vs uq for different values of the fl/uj param- 
eter. Those energy diagrams may be applied to analyse 
either magnetic resonance or dynamical localisation or 
shaken optical lattices. The zero crossing of the energy 
represent magic values where the effective magnetic en- 
ergy or quantum tunneling are frozen at values different 
from = J = 0. 



IV. RENORMALIZATION VS ujo VALUE 



Here k even, I odd, and we introduced 



A+ 



X — 4 



(17) 



Eqs. (|15p and (fTB|) are independent and a complete so- 
lution is obtained by adding the solutions of those equa- 
tions. Eqs. (jl6p may be rewritten as 



Xj-i + Xj+i 



(18) 



The energy renormalization, to be investigated on the 
basis of different theoretical approaches will be concen- 
trated on magnetic resonance case, but the analysis of 
Sec. H has demonstrated that the substitution J — ujq/2 
allows to apply our results also to the shaken optical lat- 
tices. 

a) Jo solution The two-level energy splitting at w 3> wq 
derived in Eq. ([5]) by the magnetic resonance treatment 
leads to the TZ renormalization given by j7o Bessel func- 
tion of Eq. The j7o renormalization approximation 
corresponds to the following quasienergies: 



by imposing Xj = Aj for even j, xj = Bj for odd j, with 



±^ Jo ( - 

2 \ uj 



(24) 



jAoj/n for odd j , (19a) 

for even j. (19b) 



This solution predicts a freezing for whatever ujq at the 
il/uj values corresponding to the zeros of the Jo Bessel 
function, but its validity is limited to wq ~ 0. 
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FIG. 1: (Color online) Central Brillouin zone of the 
quasienergy £+,o vs ojo, both measured in u units, for dif- 
ferent values of fl/uj, between and 5. The quasienergy e_,o 
is the opposite of e+,o- Quasienergies calculated by truncat- 
ing the continued fraction to seven terms. Freezing occurs 
when the quasienergy is equal to zero at values different from 
(jjQ — 0. For Jl/cj = 2.4 close to the Jo Bessel function first 
zero, owing to the quasienergy flatness at ojq ~ a nearly 
perfect freezing is reached in a large range of low ojq values. 
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b) S -corrected solution On the basis of a magnetic res- 
onance semiclassical treatment, refs. [13, [Hj derived an 
WQ-dependent correction to the j7o renormalization. That 
correction leads to the following quasienergies and renor- 
malization: 



± 



Jo 



S 



n 



Jo 



?ia;,(25) 
(26) 



Here S{x) is a product of Jn ordinary Bessel functions 
well approximated by the following expression (sij : 



Six) 
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2025a;'' 



[aix)j2ix) + Pix)Mx) - 7(a;) J6(a:)] 



where a{x) 
6(408- 74x2 



75 (5 - a; 
23a;Vl6) and 



2/4) x\ P{x) 



(27) 



145x2 (3 - x2/2) /49. Within the following Section 
this solution will be used for calculations around the first 
and second zeros of the zero-th order Bessel function, 
and the validity limits for the wo/w application range 
will be discussed there. 

c) Continued fraction A general approach to derive 
the 7?. renormalization coefficient is based on the Floquet 
quasienergies derived in previous Section, leading to 



hiijQ/2 



(28) 



FIG. 2: (Color online) On the top Ti vs fi/cj calculated at 
ajo/<^ = 0.1 on the basis of different approximations: (a) from 
the continued fraction containing terms up to 9 terms, the 
J'o solution of Eq. ([5]) and the 5-corrected one of Eq. (|26l) : 
(b), (c), (d) and (e) from the continued fraction limited to 
seven four, three and one terms, respectively. On the bottom 
quasienergies e+.o (continuous lines) and e„.o (dashed lines) 
vs u}o at Q./uj = 2.405. Thicker blue lines calculated on the 
basis of the four terms continued-fraction; thinner red lines 
on the S-corrected solution. The central black horizontal dot- 
dashed line based on Eq. ((24]). Diagonal lines for the ojo — > 
quasienergies. 



This equation allows a numerical determination of TZ 
without restricting to the low uq values where the JTq 
approximation (without or with the S function) is valid. 
The quasienergy can be derived cither from the continued 
fraction solution of Eq. (20), truncated to a finite num- 
ber j of terms, or from a diagonalization of the system of 
Eqs. (jl8p truncated to a finite number of equations with 
2j + 1 terms. The Hamiltonian diagonalization approach 
was applied in ref. Q for the renormalisation calculations 
in both the low-frequency and high-frequency regimes. 



V. NUMERICAL RESULTS 

The TZ renormalization coefficient is a complex func- 
tion of the system parameters loo/uj and ^/lu, and the 
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FIG. 3: (Color online) In (a) renormalization coefficient TZ 
and in (b) effective tunneling Jeff vs Q/u) ratio for different 
values of J/cj. Continuous blue lines for J/cj = 0.02, dashed 
green lines for J/ui — 0.2 and dotted red lines for J/u = 0.4. 
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FIG. 4: (Color online) [Q/u 



I freeze 



values required to freeze 



the tunneling coefficient {Jeff = 0) as a function of the J/ui 
unperturbed tunneling coefficient. Continuous lines from the 
S corrected expression and dots from the continued fraction 
solution. Closed blue dots for the first zero crossing of the 
eigenenergies and open blue ones for the second zero cross- 
ing. The stars on the horizontal axis denotes the J/lj values 
explored in the shaken optical lattice experiments. 



previous Section approaches may be used for numerical 
analyses at different parameter values. Fig 2(a) reports 
tlie R results at ujq/oj ~ 2J/uj ~ 0.1 vs fl/uj, obtained 
using Eq. ([25]) linking tliat coefficient to the Floquet 
quasienergies. For the quasienergies determined from the 
continued fraction solution, and also from the diagonal- 
ization of the system of Eqs. (fT8|). the results of Fig. 2(a) 
shows a slow convergence at low ujq values, as already 
pointed out by Autler and Townes 2J], the number of 



required terms in the continued fraction depending on 
the Q/uj value. In the wq, J — > limit TZ is well ap- 
proximated by the jQ{Q/io) function, the S correction 
vanishing there. 

For the quasienergies dependence on ujq at fl/ui = 
2.405, Fig. 2(b) compares the continued fraction solu- 
tion to the J'o solution and the S'-corrected one. The j7o 
solution leads to a horizontal line close to x-axis because 
CI/lj corresponds to the Bessel first zero, indicating that 
it approximates the quasienergies only for wo ~ 0. On 
the contrary the S'-corrected solution approximates well 
the quasienergies for a large range of parameters, at least 
for luq/uj up to 0.5 corresponding to J/cu up to 1. 

For the shaken optical lattice experiments where the 
condition wq = 2 J <g; w is satisfied, as in most cases, the 
j7o solution is well appropriate for TZ. At larger J values, 
the S function correction to Jq can be used for the full 
range of the parameters explored so far in experiments. 
Fig.[3]reports a S-correction based analysis of the shaken- 
lattice renormalization at increasing values of J/lj. For 
J/oj = 0.4 the correction to TZ shown in Fig. ^a) is ten 
percent smaller than the J/w w value, but becomes 
larger increasing J. Because the most important quan- 
tity is the tunneling coefficient itself. Fig. ^h) shows the 
Jeff « dependence on fl/uj at increasing values of J/lo. 
Notice that increasing J/oj the Jeff = freezing config- 
uration is reached at an Jl/w value lower than the Bessel 
first zero. 



Fig. H shows the [fi/a 



J freeze 



values required to pro- 



duce a tunneling freeze for a given J/uj initial value. We 
plot the values associated to the first and second zero- 
crossing of the eigenenergies, corresponding to the first 
and second zero of the Bessel function within the Jo so- 
lution. A comparison between the S'-corrected solution 
and the continued fraction solution is presented, confirm- 
ing that for most shaken-lattice experiments performed 
so far, the S'-corrected solution provides a simple and pre- 
cise determination of the modified tunneling parameter. 
For a larger range of parameters the continued fraction 
solution should be used. The data points at fi/w — >■ 
correspond to the quasienergy crossings in absence of rf 
drive and don't have a physical meaning. Notice that 
freezing can be produced also applying uj values lower 
than J, a regime was not yet examined in the experi- 
ments. It may be noticed that the general dependence of 
the freezing value of Fig. 2] is similar, although not iden- 
tical to the Bloch-Siegert shift dependence investigated 
in P, H, [13, m, [13-133 • In fact for an applied oscillating 
field, as in the present magnetic resonance configuration, 
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all the crossings and anticrossings of the energy levels 
are shifted towards lower ujq values by increasing fl [3^ . 
The Bloch-Siegert shift of the magnetic resonance is as- 
sociated to the position of the energy anticrossings, while 
the freezing point is associated to the zero crossing of the 
eigenenergies. 

VI. CONCLUSIONS 

We have examined the energy rcnormalisation of a two- 
level system usually associated to the dressed atom ap- 
proach, but also derivable from a semiclassical analysis of 
the magnetic resonance. The ^-solution derived in ref. @ 
for the Zceman freezing of optical clocks is here applied to 
the optical lattice experiments. Within that framework 
the standard zero-Bessel dependence on the amplitude 
of the electromagnetic field amplitude, valid only at zero 
magnetic field, was extended to derive a general depen- 
dence on the magnetic field amplitude. That result is 
important for the main target of the present work, to 



use the magnetic resonance results in order to perform 
an accurate analysis of the renormalisation occurring for 
the atomic quantum tunneling between the minima of 
an optical lattice in the shaken lattice experiments. The 
magnetic resonance correction to the energy renormali- 
sation allows us to derive a very general formula for the 
tunneling renormalisation in shaken optical lattices. The 
conditions for the complete cancellation of the tunnel- 
ing rate are functions of the tunneling energy without 
shaking and of the modulation frequency. The precise 
determination of the tunneling under different driving 
conditions will lead to a better control in the quantum 
simulation experiments based on optical lattices. 
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